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It is known that the dynamical evolution of a system, from an initial tensor product state of 
system and environment, to any two later times, t\,ta (t<2 > ti), are both completely positive (CP) 
but in the intermediate times between ti and £2 it need not be CP. This reveals the key to the 
Markov (if CP) and nonMarkov (if it is not CP) avataras of the intermediate dynamics. This is 
brought out here in terms of the quantum stochastic map A and the associated dynamical map B - 
without resorting to master equation approaches. We investigate these features with four examples 
which have entirely different physical origins (i) a two qubit Werner state map with time dependent 
noise parameter (ii) Phenomenological model of a recent optical experiment (Nature Physics, 7, 
931 (2011)) on the open system evolution of photon polarization, (iii) Hamiltonian dynamics of a 
qubit coupled to a bath of N qubits and (iv) two qubit unitary dynamics of Jordan et. al. (Phys. 
Rev. A 70, 052110 (2004) with initial product states of qubits. In all these models, it is shown that 
the positivity/negativity of the eigenvalues of intermediate time dynamical B map determines the 
Markov/non-Markov nature of the dynamics. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 
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I. INTRODUCTION 

Understanding the basic nature of dynamical evolu- 
tion of a quantum system, which interacts with an in- 
accessible environment, attracts growing importance in 
recent years 0, ■ This offers the key to achieve control 
over quantum systems - towards their applications in the 
emerging field of quantum computation and communi- 
cation [3(. While the overall system-environment state 
evolves unitarily, the dynamics governing the system is 
described by a completely positive (CP), trace preserving 
map [j-l6|. 

Markov approximation holds when the future dynam- 
ics depends only on the present state - and not on the 
history of the system i.e., memory effects are negligi- 
ble. The corresponding Markov dynamical map consti- 
tutes a trace preserving, CP, continuous one-parameter 
quantum semi-group 0, HI- Markov dynamics govern- 
ing the evolution of the system density matrix is con- 
ventionally described by Lindblad-Gorini-Kossakowski- 
Sudarshan (LGKS) master equation 0, Q |f = Cp, 
where C is the time- independent Lindbladian operator 
generating the underlying quantum Markov semi-group. 
Generalized Markov processes are formulated in terms of 
time-dependent Lindblad generators and the associated 
trace preserving CP dynamical map is a two-parameter 
divisible map [S| fiol ] , which too corresponds to memory- 
less Markovian evolution. 

Not completely positive (NCP) maps do make their 
presence felt in the open-system dynamics obtained from 
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the joint unitary evolution - if the system and environ- 
ment are in an initially quantum correlated state [TT14u3 j . 
In such cases, the open-system evolution turns out to be 
non-Markovian [l5[. However, the source of such non- 
Markovianity could not be attributed entirely to either 
initial system-environment correlations or their dynami- 
cal interaction or both. This issue gets refined if initial 
global state is in the tensor product form, in which case 
the sole cause of Markovianity/non-Markovianity could 
be attributed to dynamics alone. It is known that the 
time evolution of a subsystem from an initial tensor prod- 
uct form to two different later times, ti,ti {ti > ii), are 
both CP. However the dynamics in the intermediate time 
steps between t\ and ti need not be CP. The quantum 
stochastic A and dynamical B maps - first introduced as 
a quantum extension of classical stochastic dynamics - by 
Sudarshan, Mathews, Rau and Jordan (SMRJ) [g| nearly 
five decades ago, offer an elegant approach to explore 
Mar kovian/non-Markovian nature of open system evolu- 
tion. The interplay of A and B maps at intermediate 
times, to bring out the Markov or non-Markov avataras 
of open system evolution, is established in this paper. 

To place these ideas succinctly, there are three ba- 
sic aspects in open system quantum dynamics: (1) na- 
ture of dynamical interaction between the system and 
its environment, (2) role of initial correlations in system- 
environment state and (3) nature of dynamics at inter- 
mediate times. Last few years have witnessed intense 
efforts towards understanding these 0-[2l| . The third is- 
sue is the focus here to discern the Markov/non-Markov 
nature of dynamics in terms of intermediate time A and 
B maps. 

The contents are organized as follows: In Sec. II some 
basic concepts @ on A and B maps are given. The emer- 
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gence of CP/NCP maps, at intermediate times, under 
open system dynamics is discussed in Sec. III. Sec. IV is 
devoted to a powerful link (brought out by Jamilkowski 
isomorphism) between the B map and the dynamical 
state. Some illustrative examples of dynamical B map 
to investigate the CP/NCP nature of dynamics at in- 
termediate times are discussed in Sec. V. The examples 
are chosen from different origins: one based entirely from 
the general considerations of Jamiolkowski isomorphism; 
second one on the recent all-optical open system exper- 
iment to drive Markovian to non-Markovian transitions; 
the other two examples are based on open system Hamil- 
tonian dynamics. In all these four examples, no mas- 
ter equation is employed in the deduction of Markov to 
non-Markov transitions - but the CP/NCP nature of the 
intermediate dynamical map (via the sign of the eigen- 
value of the B map) has been invoked. Sec. VI has some 
concluding remarks. 

II. PRELIMINARY IDEAS ON DYNAMICAL A 
AND B MAPS 

The stochastic A and dynamical B maps @ transform 
the initial system density matrix ps(to) to final density 
matrix ps(t) via, 

[Ps(t)} blb2 = [ A (t,to)] blb2;aia2 [ps(to)] axa2 , (1) 

ai,o 2 

[ps(t)] bxb2 = [B(t,t )] biai . b2a2 [ps(t )} aia2 , (2) 

ai ,a2 

ai, a 2 , &i,i>2 = 1,2, . . .,d 
where the realligned matrix B is defined by, 

B bl a\\b'ia-i Ab 1 b 2 ;a 1 a 2 - (3) 

The requirement that the evolved density matrix ps{t) 
has unit trace and is Hermitian, positive semi-definite 
places the following conditions on A and B [(|: 

Trace preservation : ^^6 1 b 1 :« 1 n 2 —^a 1 a 2 , 

^Z B b 1 a 1 :b 1 a 2 = $ ai a 2 , 
t>l 

Hermiticity: A^j^,, = K 2 b 1 ;a 2 a 1 > 

Bb t a t ;b 2 a 2 = Bl 2a2 . biai (4) 

Positivity : ^ x* bi x b2 A blb2 . aia2 y ai y* 2 > 0, 

°1 < a 2 < b l < b 2 

X ' bl Vo.! B biai . b2 a 2 X b2 V* 2 > 

°1 < a 2 < h l < b 2 

It may be readily identified that the dynamical B map 
is positive, Hermitian d x d 2 matrix with trace d - cor- 
responding to CP evolution. We would also like to point 
out here that the composition of two stochastic A-maps, 

Ax * A 2 tranforming p s (t ) —k ps(h) — ^ Ps(h) is 
merely a matrix multiplication, whereas it is not so in 
its S-form. 

III. CP/NCP NATURE OF INTERMEDIATE 
TIME A AND B MAPS 

Let us consider unitary evolution of global system- 
environment state ps(to) ® PE(to) from an initial time 



to to a final time t 2 - passing through an intermediate 
instant t\ (i.e., to < t\ < t 2 ). The A- map associated 
with to to t\ and that between to to t 2 are identified as 
follows: 

Ty e [U(t 3 ,t )ps(t ) <E> PEitoW^t^to)] = Afc, t ) Psih) 
= p S {t j ), J =1,2. (5) 

The stochastic map A(tj,to) is completely positive (cor- 
respondingly the dynamcal B(tj,to) matrix is positive). 

In order to identify the intermediate stochastic map 
A(t2, t±), we make use of the composition law of unitary 
evolution U(t2,t ) — U(t2,ti)U(ti 7 t ): 

Tr E [t/(i 2 ,ti) {u{t u to)ps(t ) ® PB(*o)(7 t (ti,to)} U^{t 3 ,ti)] = 

A(t 2 ,t )ps(t ). (6) 

However, this does not lead naturally to A(t 2 ,t ) — 
A(t 2 , ti)A(t±, to) for the A map. Invoking Markovian 
approximation (memory-less reservoir condition [221 ]) 

{U(t 1 ,t )ps(t ) ® PEitoW^t^to)} « p S (h) ® pjsfa), 

the LHS of ([5]) may be expressed as, 

Tr £ [U(h,ti)ps(h) <8> p E {ti)U\t 2 , h)] = A(t 2 ,h) p s (h). 

(7) 

Further, substituting j = 1 in ([5]1 and expressing 
PS (to) = A~ 1 (ti, t )ps(ti) in ^ the intermediate A map 
A(t 2 ,t L ) is identified: 

A(t 2 ,t 1 ) = A(t 2 ,to) A ~ 1 (t l ,to). (8) 

In other words, when the environment is passive ( Marko- 
vian dynamics), the intermediate A- map has the divisi- 
ble composition as in flSJ). In such cases A(ti,t 2 ) is en- 
sured to be CP - otherwise it is NCP, and hence non- 
Markovian. Correspondingly, the intermediate B-map 
B(t 2 ,ti) is positive if the dynamics is Markovian; nega- 
tive eigenvalues of B(t 2 ,ti) imply non-Mar kovianity. 



IV. THE B MAP AND THE JAMIOLKOWSKI 
ISOMORPHISM 

The Jamiolkowski isomorphism Q provides an insight 
that the £?-map is directly related to a d 2 x d 2 system- 
ancilla bipartite density matrix. More specifically, the 
action of the map A ld ® A on the maximally entangled 
system-ancilla state |V>me) = Yli=o N) *) results in the 
density matrix p ab which may be identified to be related 
to the dynamical B map i.e., 

Pab = [A M ® A] |Vme)(V>me| -> 2 B (9) 

gives an explicit matrix representation for the B-map 
(Here A Id is the identity A-map, which leaves the ancilla 
undisturbed). 

In detail, we have, 
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or (Pba)h iai .i J2a2 — d Bb iai -b 2a2 . 

In other words, Jamiolkowski isomorphism maps ev- 
ery completely positive dynamical map B acting on d di- 
mensional space to a positive definite d 2 x d 2 bipartite 
density matrix p ab (See Eq. (jlOp - ) - whose partial trace 
(over the first subsystem - as seen from the trace preser- 
vation property on dynamical map B (as in Eq.Q) is 
a maximally disordered state. One such set of bipar- 
tite d x d density matrices belong to the class that are 
invariant under U ®U [23[ - which constitute the well- 
known Werner density matrices. One may now identify 
several toy models of dynamical B maps - including the 
two qubit Werner state example motivated by the above 
remark - to investigate the nature of intermediate time 
dynamics. 

In the next section we present specific examples cho- 
sen to illustrate the features of intermediate dynamical 
maps: (i) A toy model map inspired by Jamiolkowski 
isomorphism - which is not based on any Hamiltonian 
underpinning, (ii) Recent optical experiment by Liu et. 
al., [24j on open system evolution of photon polarization 
to bring out non-Markovianity features is reinterpreted 
in terms of NCP nature of the intermediate dynamical 
map. (iii) Intermediate dynamical map in the Hamilto- 
nian evolution of a two-level system coupled to N two- 
level systems [2l| (iv) open system dynamics arising from 
a two qubit unitary evolution [ill ]. 



V. EXAMPLES 



A. A toy model dynamical map 



ied. Its use here as a valid i?-map is novel in identifying 
transitions between Markovianity and non-Markovianity 
in the dynamics as captured from their intermediate time 
behavior. 

On evaluating the corresponding A map A(t, 0) (ex- 
pressed in the standard {|0, 0), |0, 1), |1, 0), |1, 1)} basis) 



A(t,0) = 
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one can obtain the intermediate dynamical map 
A(t 2 , h) = A(t 2 ., 0)/l _1 (ti, 0). The intermediate time 
B-map B(t 2 ,ti) is given by 



B(t 2 



' 2 



p{h) 



2p(t 2 



(12) 



and 



Its eigenvalues are Ai = A 2 = A3 = \ \ \ — 

x _ 1 {1 , 3p(t 2 ) \ 

A choice p(t) = cos 2M (at) for any M > 1 leads to 
NCPness of the intermdiate map - as the eigenvalues 
A1.2.3 = A of B(t 2 ,t\) may assume negative values - 
and hence non-Markovian dynamics ensues. We have 
plotted the negative eigenvalue A of B(t 2 - ti) as a function 
of n = t-xjt\ and for typical values of M = 1,3,5 in 
Fig. 1. This reveals transitions from Markovianity to 
non-Markovianity and back in this model. 



The two qubit Werner density matrix is a natural 
choice for a prototype of dynamical i?-map - arising 
from general considerations of the Jamiolkowoski isomor- 
phism: 



B(t,0) 



[1 -*>(*)] 



2 1 



(11) 



with a time dependent noise parameter < p(t) < 1, 
and = -± (|0,0) - |1,1)) is the Bell state. For a 

dynamical map, time dependence in p(t) occurs due to 
the underlying Hamiltonian evolution. This state is espe- 
cially important in that it exhibits both separable and en- 
tangled states, as its characteristic parameter p(t) is var- 
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FIG. 1: A plot of the eigenvalue A of B(t2, ti) versus /j, — I 2 - 
for different values of M. The dynamics is non-Markovian 
when A assumes negative values and otherwise it is Marko- 
vian. 
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Another choice pit) = e~ at corresponds to a CP in- 
termediate map - resulting entirely in a Markovian pro- 
cess. In this case, we also find that Afa, t\) = Afo — tx) 
and this forms a Markov semigroup. However, if p(t) — 

e~ at , (j3 7^ 1), the intermediate map is still CP (and 
hence Markovian), though A(t2,ti) ^ A(t2 — t\) and 
therefore, it does not constitute a one-parameter semi- 
group. 

Furthermore, we wish to illustrate through this toy 
model that concurrence of p a b(t) = ^£?(i,0) (given by 
C = 3p ^~ 1 ) can never increase as a result of Markovian 
evolution. This is because ensuing dynamics is a local CP 
map on the system. Any temporary regain of system- 
ancilla entanglement during the course of evolution is 
clearly attributed to the back-flow from environment to 
the system - which is a signature of non-Markovian pro- 
cess. This feature is displayed in Fig. 2 by plotting the 
concurrence of p a b(t) for different choices of p(t). 




at 



FIG. 2: Concurrence C = 



3p(t)- 



of the system-ancilla state 



p a b(t) = [1 ~4 (t)1 h ® h +p(t) |V'me)(V'me|, vs scaled time at, 
for the following choices (i) Markov process: p(t) = e~ at (solid 
line) and (ii) non-Markov process: p(t) = cos 2M (at), M = 1 
(dashed line) and M = 5 (dot-dashed line). Note that there is 
a death and re-birth of entanglement (dash, dot-dashed lines) 
due of back-flow from environment. 



B. Optical Experiment 

Recently, Liu et al [24j reported an optical experiment 
on the open quantum system constituted by the polar- 
ization degree of freedom of photons (system) coupled to 
the frequency degree of freedom (environment) . They re- 
ported transition between Markovian and non-Markovian 
regimes. 

The dynamical evolution of the horizontal and vertical 
poloarization states \H), \ V) of the photon is captured by 
the following transformation: 



\H)(H) h+ \H)(H\ 

\V)(V) h+ \V)(V\ 

\H)(V) h+ K*(t)\H)(V\ 

\V)(H) ^ n(t)\V)(H\ 



(13) 



where denotes the decoherence function, magnitude 
of which is modelled as (for details see [24| ) , 



|«(t)| 



l-AAi(l-Ax) sin 2 (r Aw) (14) 
< A x < 1. 



The corresponding A and B maps (in the 
{HH, HV, VH, VV} basis) are readility identified 
to be, 



A(t,0) 



B(t,0) = 



( 1 

K*(t) 

n(t) 

Vo 1 

( 1 K*(t) 



V n(t) 1 



(15) 



We construct the intermediate time dynamical 
map B(t2,tx) from the corresponding A(t2,tx) = 
A(t 2 ,0)A~ 1 {t 1 ,0) to obtain, 



B(t 2 ,tx 



/ 1 


























k(* 2 ) 









K*(*i) 




1 



*2) \ 
t, 1 \ 



(16) 



Eigenvalues of B(t2, t\) are given by, 

«(*2) 



Ai. 



1± 



K(tl) 



2,3 



0. 



(17) 



The eigenvalue A4 can assume negative values indicating 
Markovian/non-Markovian regimes. A plot of the neg- 
ative eigenvalue as a function of Ax , for different ratios 
tz/tx, is given in Fig. 3 - where one can clearly see the 
Markovian (A4 > 0) and non-Markovian (A4 < 0) regimes 
in this model. 




FIG. 3: A plot of the eigenvalue A4 of B(t2, ti) versus Ai for 
different values of fi — |^-. 



C. Hamiltonian evolution of a two level system 
coupled to a bath of N spins 

We now present a Hamiltonian model, which give rise 
to explicit structure of time dependence in the open sys- 
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tern evolution. Interaction Hamiltonian considered here 

is m 



H = 



A 



N 



Ok; 



(18) 



This is a simplified model of a hyperfine interaction of 
a spin- 1/2 system with N spin- 1/2 nuclear environment 
in a quantum dot. Taking the initial system-environment 
state to be ps(0)<S> the dynamical ^4-map is obtained 



by evaluating Tr B U(t,0)p s (0)< 
U(t,0) = Exp[-iHt]): 



W(t,0) (where 



1 1 

A(t, 0) = -{\-x{t))a z ®a z + -{l + x{t))I 2 ®I 2l 



x(t) = cos iV —= 



2At 



(19) 



From this, the intermediate map A(t2,ti) (see ©) and 
in turn the corresponding B(t2,ti) may be readily eval- 
uated. We obtain, 



B(t 2 , h) 



x(t 2 



I? 



{u 3 



iv y ). (20) 



The eigenvalues of B (t 2 , 1 i ) are0,0,l± . Clearly, the 
intermediate time dynamics exhibits NCP as one of the 
eigenvalues i.e., A = 1 — °f B(t 2 , t\) can assume neg- 
ative values. We illustrate regimes of Markovianity/non- 
Markovianity revealed via positive/negative values of A 
(plotted as a function of p = tijti) in Fig. 4. 
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on the system- environment initial state pse(0) — Ps(0)<8) 
p E {0) = | (J 2 + a x ) ® i (I 2 + a z ) . The A(t, 0) map is 
given by, 



A{t,ti) = ~ (l + cos(wi)) I 2 ®I 2 +~ (1 - cos(wt)) 0- z <g)cr*. 



Following (O, we obtain 



(22) 



= ^ ( J 2 ® J 2 + CT 2 



cos(wi 2 ) , 
2 cos(wii) 



a„) . (23) 



The eigenvalues of the B-map are given by 0, 0, 1 ± 



COS U)t2 
COS Ujt\ 



The eigenvalue A = 1 ± 



COS (jjt2 



can assume neg- 

cosLjri ° 

ative values - bringing out the non-Markovian features 
prevalent in the dynamical process. Fig. 5 illustrates the 
transitions from Markovianity to non-Markovianity. This 
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FIG. 5: The plot of the eigenvalue A = 1 - as a 

r ° cos (jjti 

function of p = ^ ■ The periodic transitions of A from positive 
to negative values indicates the transition of the process from 
Markovian to non-Markovian. 

model, with initially correlated states, has been explored 
before in Refs. [HI [H[ and the dynamical map turned 
out to be NCP throughout not merely in the intermedi- 
ate time interval). 



FIG. 4: The variation of the eigenvalue A of B(t2, ti) (as a 
function of p — fe/ti) from positive to negative values and 
back with the passage of time for different values of N. 



D. Two qubit unitary evolution 

We now consider the open system dynamics arising 
from the unitary evolution [ll| 



U(t, 0) 



cos(w t/2) I 2 ® J 2 — i sin(w t/2)a z 



(21) 



VI. SUMMARY 

In conclusion, a few remarks on a variety of defini- 
tions of non-Markovianity in the recent literature may 
be recalled here. Mainly the focus has been towards cap- 
turing the violation of semi-group property [TBI . [l7j or 
more recently - its two-parameter generalization viz the 
divisibility of the dynamical map [9j. llQ|. Yet another 
measure, where non-Markovianity [l8[ is attributed to 
increase of distinguishability of any pairs of states (as a 
result of the partial back-flow of information from the 
environment into the system) and is quantified in terms 
of trace distance of the states. It has been shown that 
the two different measures of non-Markovianity - one 
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based on the divisibility of the dynamical map [To| and 
the other based upon the distinguishability of quan tum 
states [l8| - need not agree with each other |l9j . A 
modified version of the criterion of Ref. fTp| was pro- 
posed recently 21]. In this paper we have established 
the interplay of stochastic A and dynamical B maps at 
intermediate times, revealing Markovian/non-Markovian 
regimes. We have explored four different examples re- 
vealing the features of intermediate time maps originat- 
ing from variety of physical mechanisms : (i) A toy 
model map inspired by general considerations based on 
Jamiolkowski isomorphism - which explores a two qubit 
Werner state with time-dependent noise parameter as a 
dynamical map (ii) A reinterpretation of the phenomeno- 
logical model explaining the recent optical experiment 
by Liu et. al., 24] in terms of NCP nature of the in- 
termediate B map. (iii) Hamiltonian evolution describ- 



ing the hyperfine interaction of a spin-1/2 system with 
N spin-1/2 nuclear environment in a quantum dot [2lj ] 
displaying Markovian/non-Markovian behaviour and (iv) 
Unitary evolution of Jordan et. al., [TlT | - wherein ini- 
tial system-environment two qubit is chosen in a product 
state. Here too, intermediate time dynamical map ex- 
hibits Markov/non-Markov regimes. It is interesting to 
note that the dynamics had been identified to be NCP 
throughout not merely in the intermediate time interval 
- when initially correlated states were employed fill . [l5[ . 
Placing these two results together, brings forth that the 
source of non-Markovianity in this model is attributable 
entirely to the unitary dynamics — rather than initial 
correlations of system-environment qubits. We have thus 
exposed the underlying features of intermediate time A 
and B maps to bring out clearly if the dynamics relies on 
past history of the states or not. 
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